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Abstract 

A  stable  Luenberger  Observer  was  developed  to  reconstruct  the 
error  states  in  modal  coordinates  for  a  satellite  in  a  oeriodic 
orb  it  about  the  Earth-Moon  Lagrangian  point,  L3.  Observer  Dole 
placement  for  the  third  root  was  found  to  be  linear  with  applied 
gain.  Convergence  of  the  observer  estimate  to  the  true  state 
vector  was  rapid  in  all  cases,  but  depended  upon  the  initial  con¬ 
ditions  considered.  The  addition  of  the  observer  did  not  impede 
the  long-term  behavior  of  the  controller  in  stabilizing  the  orbit. 
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CLOSED-LOOP  CONTROL  OF  A  SATELLITE 


IN  AN 

UNSTABLE  PERIODIC  ORBIT 
ABOUT  L3 


I.  Introduction 


Background 

The  n-body  problem  came  into  existence  when  celestial  mechanicians 
began  studying  the  solar  system  with  an  eve  toward  evaluating  the 
equations  of  motion  of  the  various  gravitationallv-attracted  bodies 
found  therein.  The  motion  of  two  bodies  in  a  circular  orbit  about 
each  other  and  affecting  each  other  only  through  their  mutual  forces 
of  gravitation  is  now  readily  understood,  and  when  a  third  bodv, 
infinitesimal  and  massless  in  comparison  to  the  others  is  added  to 
the  system,  the  restricted  three-bodv  problem  is  established.  The 
discovery  of  the  Trojan  group  of  satellites  in  the  Sun-Jupiter  system 
(1906)  confirmed  Comte  Joseph  Louis  Lagrange's  1772  mathematical 
depiction  of  such  a  configuration.  His  work  developed  the  existence 
of  equilibrium  solutions  to  the  equations  of  motion  where  the  gravi¬ 
tational  forces  on  the  infinitesimal  third  bodv  as  exerted  bv  the 
other  two  would  exactly  balance,  thereby  rendering  it  motionless 
with  respect  to  a  rotating  coordinate  reference  frame.  These  five 
equilibrium  positions,  LI  through  L5,  called  alternatelv  Lagrangian 
or  libration  points,  also  exist  for  the  Earth-Moon  svstem  and  are 
shown  in  Figure  1.  For  a  rotating  coordinate  frame  of  reference 
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centered  at  the  Earth-Moon  center  of  mass,  the  Lagrangian  points 
are  five  fixed  points  at  which  a  third  mass  would  be  stationary  if 
placed  with  zero  velocity  relative  to  the  frame.  Through  infinites¬ 
imal  analysis  (Ref  17:427)  it  can  be  shown  that  motion  in  the 
neighborhood  of  the  collinear  points  (LI,  L2,  and  L3)  is  inherently 
unstable  since  a  positive  real  eigenvalue  exists  in  each  case. 

Similar  evaluation  yields  motion  in  the  vicinity  of  L4  and  L5  to 
be  infinitesimally  stable.  In  1907,  Charlier  (Ref  22)  analyzed  two 
classes  of  periodic  orbits  around  these  triangular  points,  long  and 
short,  depending  on  the  period  or  mass  ratio  of  the  two  large  masses 
The  desirability  of  placing  satellites  at  any  or  all  of  these 
Lagrangian  points  coupled  with  the  related  instabilities  has  justi¬ 
fied  the  vast  amount  of  research  which  has  surfaced  in  the  past  two 
decades  dealing  with  the  orbital  control  of  this  system. 

The  major  asset  attained  by  satellites  placed  at  the  Lagrangian 
points  is  that  they  remain  in  a  motionless  configuration  with  respect 
to  the  Earth  and  the  Moon  at  all  times.  Many  uses  for  which  such 
positioning  would  be  ideal  have  been  proposed.  One  is  in  the  develop 
ment  of  a  communications  link  between  the  Earth  and  semipermanent 
bases  on  the  far  side  of  the  Moon.  Relay  satellites  placed  at  LI  and 
L2  could  provide  uninterrupted  communications  to  the  Earth  as  well 
as  to  far-reaching  lunar  exploration  vehicles.  Communication  with 
interplanetary  probes  would  be  similarly  facilitated  since  the  Earth 
as  a  noise  source  could  thus  be  virtually  eliminated.  Instead  of 
a  lunar  parking  orbit,  a  station  at  L2  would  provide  an  infinite 
launch  window  for  lunar  activities  since  it  is  stationary  with  respec 
to  the  surface.  Telescopes  of  all  types  would  be  free  from  the 


restricting  layers  of  the  Earth's  atmosphere,  or,  depending  on 
position,  be  entirely  blocked  from  Earth  interference  by  the  Moon. 
Situated  outside  of  the  magnetosphere,  satellites  to  study  solar 
plasma,  cosmic  rays,  and  other  cosmic  phenomena  would  find  ideal 
settings,  similar  to  the  ISEE-3  research  satellite  already  in  use. 

L4  or  L5  satellites  could  provide  long-term  solar  flare  observations, 
early  warning  programs,  and  secure  the  previously  unaccomplished  task 
of  performing  simultaneous  measurements  at  different  locations 
within  cislunar  space,  thus  separating  temporal  from  spatial  varia¬ 
tions  for  the  first  time  (Refs  5,  6,  8).  Militarily,  a  defense 
satellite  network  placed  in  orbit  around  L3,  L4,  and  L5  would  pro¬ 
vide  continual  global  coverage,  especially  vital  in  today's  potentially 
explosive  environment.  Synchronous  satellites  are  not  totally  efficient 
in  such  instances  since  they  exist  in  equatorial  orbits,  leaving  the 
polar  regions  essentially  hidden.  For  the  Lagrangian  system,  this 
restriction  does  not  apply,  since  these  points  lie  in  the  Earth-Moon 
plane  rather  than  equatorial  (Ref  26).  Also,  they  are  distant  enough 
that  the  lead  time  for  any  interceptor  is  sufficient  for  adequate 
evasive  maneuvers  to  be  implemented.  Prior  work  by  Captain  Tilton 
(Ref  23)  established  the  existence  of  appropriate  stable  orbits  about 
L4  and  L5  in  which  reconnaissance  satellites  could  be  placed  with 
relatively  little  cost  and  remain  for  up  to  fifty  years.  Only  the 
lack  of  a  similarly  stable  orbit  about  L3  prevents  the  completion  of 
the  triad  for  total  coverage  and  is  the  subject  of  this  and  other 
research  endeavors. 
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Svstem 


The  restricted  four-body  problem  is  a  more  realistic  extension 
of  the  three-body  problem  for  evaluation  since  the  large  perturbative 
effects  of  the  Sun  are  no  longer  neglected.  However,  in  this  case, 
the  Hamiltonian  becomes  an  explicit  function  of  time  and  there  is  no 
longer  a  constant  of  the  motion.  It  has  been  shown,  however,  that 
families  of  periodic  orbits  still  exist  about  the  five  T.igrangian 
points . 

The  modelling  assumptions  used  in  previous  research  were  maintained 
here.  The  Earth  is  taken  to  be  moving  in  a  circular  orbit  about  the 
Sun  with  a  period  of  365.256365  days.  The  Moon's  motion  is  approximated 
as  a  periodic  orbit  about  the  Earth  as  developed  by  Professor  Wiesel 
in  his  work  (Ref  25),  which  is  more  realistic  than  a  circular  one. 

All  bodies  are  considered  to  lie  in  a  plane,  so  that  perturbations  due 
to  inclination  and  eccentricity  effects  have  been  neglected.  The 
importance  of  these  are  described  in  Reference  4.  In  later  cases, 
where  greater  stability  characteristics  were  sought  for  the  nonplanar 
(z)  mode,  inclination  and  eccentricity  had  to  be  included.  Effects 
due  to  the  presence  of  other  solar  system  bodies  and  the  Sun's  radia¬ 
tive  pressure  may  be  considered  negligible  and  have  been  disregarded. 

Coordinate  Frame 

Even  though  there  is  no  longer  a  constant  of  the  motion,  which 
instigated  the  use  of  a  rotating  reference  frame  with  respect  to  the 
restricted  three-body  problem,  that  configuration  is  maintained  here. 

The  rectangular  reference  frame,  depicted  in  Figure  2,  has  its  origin 
at  the  Earth-Moon  center  of  mass  and  one  axis  coinciding  with  the  line 
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connecting  the  Earth  and  the  Moon.  This  system  rotates  at  the  average 
angular  velocity  of  the  Moon  about  the  Earth,  with  the  Lagrangian 
points  maintaining  their  relative  positions  within  the  fram>.  at  all 
times. 


Periodic  Orbit 

With  the  satellite  assumed  massless  for  the  sake  of  calculation 
and  the  attracting  bodies  having  their  respective  masses,  the  Hamil¬ 
tonian  for  the  system  was  derived  by  Captain  Shelton  to  be 


H  =  %(P  2  +  P  2+  P  2)  +  <3(P  R  -PR) 
x  y  z  x  y  y  x 

•  • 

+  g3  _wo  )p  +  (rp  +  5>P  )P 


_  _  Glml 


\^~r  +r|  |r  -  | 

mi+Tn2  m  mi+m2  m 


where  Mo=m1+m2+m3  is  the  sum  of  the  masses  of  the  Earth,  Moon,  and 

Sun  respectively,  the  scale  factor  locates  the  inertial  point  of 

o 

the  system,  ffi  is  the  angular  velocity  of  the  rotating  frame, 

R,  p,  and  r  are  the  position  vectors  from  the  origin  to  the  satellite, 
m 

Sun,  and  Moon  respectively,  with  subscripts  x,  y,  and  z  denoting  the 
components  in  those  directions,  the  superscript  r  represents  differen¬ 
tiation  in  the  rotating  frame,  Gj  is  the  gravitational  constant,  and 

P  ,  P  ,  and  P  denote  the  momenta  of  the  Hamiltonian  (Ref  20:14). 
x  y  z 

Using  the  equations  of  motion  as  derived  from  the  Hamiltonian 
and  the  theory  described  in  Chapter  II  of  Captain  Shelton's  work,  he 
was  able  to  construct  a  viable  periodic  orbit  about  L3.  This  orbit  is 
shown  in  Figure  3  and  the  initial  conditions  incorporated  in  its 
implementation  are  given  in  Table  I.  Since  such  an  orbit  was  previously 
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shown  to  be  inherently  unstable,  some  form  of  control  law  is  necessarv 
if  a  satellite  is  to  remain  in  the  vicinity  of  L3  for  any  considerable 
length  of  time. 

State  variable  feedback  control  designs  rely  on  the  availabilitv  of 
the  system's  state  vector  through  direct  measurement  or  derivation 
from  the  output.  In  this  case,  where  only  two  of  the  six  state  vector 
elements  are  accessible  through  range  and  range- rate  measurements,  a 
Luenberger  Observer  is  introduced  to  reconstruct  the  missing  states. 

This  approximate  state  vector  can  then  be  substituted  into  the  control 
law.  Rapid  convergence  to  the  true  states  is  assured  by  placing  the 
observer  poles  to  the  left  of  those  of  the  system.  Conceivablv, 
approximately  instantaneous  convergence  could  be  attained  if  the  poles 
were  specified  at  negative  infinity,  but  this  tends  to  make  the  observer 
behave  like  a  differentiator,  (Ref  15),  greatly  influenced  by  noise, 
an  undesirable  characteristic.  These  two  constraints  limit  effective 
observer  pole  placement. 


Table  I.  L3  Orbit  Initial  Conditions  (Ref  20:20) 


x  position  (km) 
y  position  (km) 
z  position  (km) 
x  momentum  (km/sec) 
y  momentum  (km/ sec) 
z  momentum  (km/sec) 


-383,537.27 

0.0 

0.0 

0.0 

-30.825537 

0.0 


II.  Floguet  Theory 


The  equations  of  motion  for  the  restricted  four-body  problem  as 
derived  from  the  Hamiltonian  (1)  are  time-varying  and  exhibit  some 
periodic  solutions.  This  means  that  the  motion  repeats  itself  pre¬ 
dictably  so  that  only  one  orbital  period  need  be  examined  to  determine 
characteristic  behavior.  In  general  form  the  equations  of  motion  may 
be  written: 

x(t)  =  f[x(t),t]  (2) 

where  x(t)  is  the  state  vector  composed  of  the  coordinates  and  momenta 
of  the  system.  Allowing  6x(t)  to  be  a  small  displacement  from  the  pre¬ 
viously  determined  periodic  orbit  taken  as  a  reference, 

<5x(t)  =  x ( t )  -  ^(t)  (3) 

the  state  vector  considered  throughout  this  study  actually  becomes 


x(t)  =  jj5x(t)  <Spx(t)  5y(t)  <5Py(t) 
Substituting  (3)  into  (2),  expanding 
reference,  and  then  linearizing  (Ref 
of  motion  written  as: 

6x(t)  =  A(t)6x(t) 


6z(t)  «Pz(t)Jr  (4) 

in  a  Taylor  series  about  the 
20:15-17)  results  in  the  equations 

(5) 


where  A(t) 


3f 

3x(t) 


is  termed  the  system  plant  matrix.  The 

x  (t) 

P 


elements  of  A(t) ,  A_(t),  are  continuous  periodic  functions  of  time 
with  period  x  =  29.530589  days,  so  that 

A(t  +  x)  =  A(t)  (6) 

defines  the  periodicity.  The  n  linearly  independent  solutions  to 
(5),  5x^(0  (i  =  1,2, ...n),  may  be  collected  to  form  an  n  x  n  matrix, 
4*  ( t )  =  |j5xi(t)  !  6x2  (t)  !  ...  !  6xn  C  t  (7) 

called  a  fundamental  matrix. 


10 


11 


<f(t  +  x)  =  M4>(  t)  (15) 

From  this  equation  at  t*0  and  recalling  (11) ,  it  is  seen  that  the 
monodromy  matrix  is  actually  equivalent  to  the  state  transition  matrix 
from  t»0  to  t=x. 

The  eigenvalues  of  M  (or  $(t,0)),  (i=l,2, . . .n) ,  are  termed 

the  characteristic  multipliers  of  A(t).  There  also  exists  a  set  of 
A  (i=l,2,...n),  related  to  the  characteristic  multipliers  by 

Xi  =  \  ln(ai) 

These  are  called  Poincare  exponents  and  may  have  real  or  complex 
values.  Since  previous  work  established  a  multiplicity  of  1  for  the 
X  they  may  be  arranged  in  the  following  Jordan  canonical  form 
(Ref  17:267):^ 

X 1  0 

A' 


0  n 

For  complex  conjugate  values,  this  format  becomes  unwieldy  since 
imaginary  vectors  would  then  be  involved  in  computations.  To  avoid  this 
inconvenience,  the  form  shown  in  (16)  may  be  used  instead.  Here,  each 
complex  conjugate  pair  forms  an  off-diagonal  2x2  submatrix  within  the 
basically  diagonal  structure. 

aj  -u>i 
u)i  <?! 


J  = 


a2 

0)2 


a? 


a  -u) 

a  n 


(16) 


where  (i  =  1,2, . 


,n)  are  the  real  components  and  ^  are  the  imaginary 


components  of  the  A^  . 
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For  purposes  of  stability,  all  Poincare  exponents  must  lie  on 
the  negative  real  axis  of  the  s-plane,  or,  if  complex,  have  a  negative 
real  part.  The  imaginary  component  yields  oscillation  in  that  particular 
related  state. 

The  PoincarS  exponents  for  the  uncontrolled  periodic  orbit  about 
L3  are  listed  in  Table  II.  These  determine  three  modes  of  oscillation, 
two  planar  and  one  nonplanar,  of  which  one  is  real  and  the  other  two 
consist  of  complex  conjugate  pairs.  As  can  be  seen,  the  third  value 
has  a  positive  real  part,  thus  rendering  the  entire  system  unstable. 

This,  then,  is  the  element  which  must  be  controlled  in  order  to  make 
this  a  useful  and  predictable  system. 

To  determine  the  actual  solution  to  (12) ,  the  fundamental  matrix 
may  be  expressed  as 

*(t)  =  0(t)eJt  (17) 

where  Q(t)  is  a  periodic  matrix  satisfying 

Q(t  +  t)  =  0(t) 

and  J  is  the  Jordan  canonical  form  associated  with  A(t)  as  defined 
above.  Substituting  (17)  into  (12), 

d  jq(t)eJtl  =  A(t)  Q(t)eJt 

dtL  -> 

differentiating, 

Q(t)eJt  +  0(t)JeJt  =  A(t)  0(t)eJt 

and  cancelling  the  e^C  terms  since  thev  can  never  be  zero,  produces 

Table  II.  Poincar4  Exponents  for  the  Uncontrolled  System 

Planar  Mode  (x)  0.0  -  1.0392i 

0.0  +  1.0392i 

Planar  Mode  ( v)  2.3921  +  O.Oi 

-2.3921  +  O.Oi 
0.0  -  1 . 1 247 i 
0.0+1 . 1247i 


Nonplanar  Mode  (z) 


Q(t)  +  0(t)J  =  A(t)0(t) 

Defining  Aft)  as  the  0(t)  matrix  partitioned  into  columns 

A(t)  =  [oj(t)  i  02(t)  On(t)J  (18) 

results  in  a  simple,  analytical  equation  for  J: 

J  =  A_1(t)[A(t)A(t)  -  Aft)]  (19) 

Now,  substituting  $(t,t  )  into  (12)  gives  the  result 

<Kt,t  )  =  A(t)$(t,t  ) 
o  o 

This  may  be  integrated  numerically  with  the  initial  conditions 
described  by  (9)  to  yield  information  regarding  a  particular  solution's 
local  variation  with  initial  conditions. 

Picking  6x^(t)  as  an  arbitrary  solution  to  (5)  and  extracting 
this  from  (17)  results  in 

6x±(t)  =  q.(t)eXiC  (20) 

From  (8)  at  t  =  0,  it  is  seen  that 

<5xi(t)  =  9(t,0)5x.(0)  (21) 

Equating  (20)  and  (21), 

5xi(t)  *  Q1(t)eXit  =  <J(t,0)  fx^(O) 
and  substituting  Q^(0)  for  6x^(0)  from  (20)  evaluated  at  t=0 
gives 

6^(0  =  Qi(t)e?'it  =  $(t,0)01(0)  (22) 

Because  0(t)  is  periodic,  0^(t)  =  0^(0), 

5x.(i)  =  O^t^V  -  O.(0)eXiT  =  <Kt,0)0.(0)  (23) 

when  (20)  and  (22)  are  equated  at  the  end  of  one  period. 

Incorporating  M  =  $(1,0)  in  the  latter  part  of  (23)  and  simpli¬ 


fying  yields  the  eigenvalue/eigenvector  problem, 


where  e^’  are  the  eigenvalues  and  the  0^(0)  the  respective  eigen¬ 
vectors  of  the  monodromy  matrix. 

For  the  complete  solution,  the  Ch(t)  must  be  computed  over  an 
entire  period.  Equation  (19)  is  essentially  an  ordinary  differential 
equation  for  A(t),  the  partitioned  0(t)  matrix,  having  as  initial 
conditions  the  eigenvectors  of  M. 
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III.  Modal  Control 


Due  to  the  system's  unstable  characteristics,  some  form  of  control 
is  necessary.  Rate  and  position  feedback  laws  were  unsuccessful  since 
all  modes,  the  stable  as  well  as  the  unstable,  were  excited  (Ref  21) 
with  appreciable  control  costs.  In  order  to  avoid  this,  Captain  Shelton 
implemented  a  modal  controller  which  resulted  in  moving  only  the  third 
root,  leaving  the  other  modes  and  the  conjugate  root  of  the  unstable 
mode  unchanged.  Although  this  method  is  not  the  only  alternative 
available,  it  is  commensurate  with  present  knowledge  in  the  control  of 
periodic  systems. 

So  that  a  controller  may  be  added  to  the  system,  the  state  equations, 
(5) ,  must  be  augmented  with  a  control  term 

B(t)u(t)  (25) 

where  u(t)  is  the  control  vector  and  B(t)  is  the  distribution  matrix 
which  determines  to  which  states  control  is  applied. 

Sx(t)  =  A(t)6x(t)  +  B(t)u(t)  (26) 

For  Captain  Shelton's  work  B(t)  was  taken  as  [o  1  0  1  o  o]T 
signifying  changes  in  momenta  in  the  x  and  y  directions  due  to  thrust 
application . 

So  that  only  the  desired  mode  is  affected  by  the  controller,  the 
system  must  be  converted  to  a  set  of  diagonal  equations.  Implementing 
the  transformation, 

5x(t)  =  A(t)n(t)  (27) 

where  n(t)  is  the  modal  vector  and  A(t)  is  as  defined  in  (18),  and 
its  derivative, 

<5x(t)  =  A(t)n(t)  +  A(t)n(t)  (28) 

the  modal  state  equations  become 

A(t)n(t)  +  A(t)n(t)  =  A(t)A(t)n(c)  +  B(t)u(t)  (2°) 
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Rearranging  (29)  produces 

n(t)  =  A-1(t)[A(t)  ‘  (t)  -  A(t)]n(t)  +  A_1(t)B(t)u(t)  HO) 
Recalling  (19) 

J  =  A_1(t)[A(c)A(t)  -  A ( c )]  (19) 

where  J  is  the  diagonal  matrix  of  Poincare  exponents  defined  pre¬ 
viously,  the  differential  equations  of  motion  in  modal  variables  are: 

n(t)  =  Jn(t)  +  A-' (t)B(t)uft)  (31) 

which  can  be  expanded  into  the  six  equations  representing  the  modal 
system. 

Previous  work  maintained  a  scalar  control  law  for  simplicitv.  So 
that  the  only  unstable  component  of  the  modal  vector  (n3)  would  be 
influenced,  the  control  law  was  defined  as 

u  =  cn3(t)  (32) 

where  G  is  the  constant  gain  term.  Transforming  this  back  to  phvsical 
variables  through  the  inverse  of  (27)  gives 

u  =  GA3_1 ( t) 6x( t)  (33) 

where  A3  *(t)  is  defined  as  the  third  row  of  the  J(t)  matrix. 

Expanding  (31)  fully  and  substituting  for  u(t)  results  in  a  linear, 
first-order,  time-varying  differential  equation  for  n  3  ( t ) : 

n  3  ( t)  =£'3  +  (A3i;(t)  +  A  3  4  ( t ))  n  3  ( t )  (34) 

where  the  subscripts  refer  to  the  rows  and  columns  resnectivelv  of 

-1 

A  (t).  Solving  this  equation  with  an  integrating  factor  and  separating 

"I  “1 

A32(t)  and  A34(t)  into  their  respective  Fourier  series  for  analvsis 
resulted  in  constant  and  periodic  components  of  each  term.  This  gave 
rise  to  a  linear  relationship  between  the  actual  and  desired  pole 
locations, 
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(35) 


"  X3  +  GYo 

where  X3  is  the  present  third  Poincare  exponent,  X3  is  its  desired 
position,  and  yq  is  the  constant  component  yielded  bv  the  Fourier 
analysis  (See  Ref  20:  31-33,  for  detailed  derivation).  Thus, 

this  procedure  reduces  to  a  pole  placement  technique. 

By  varying  the  gain.  Captain  Shelton  was  able  to  determine  that 
a  value  of  G=  .259  was  -just  adequate  to  bring  the  third  Poincare 
exponent  into  the  left-half  plane,  therbv  resulting  in  a  controlled 
stable  system.  Other  gains  simply  vary  its  position  linearlv  in 
accordance  with  the  above  law. 


1« 


IV.  Observer  Theory 

The  observability  of  a  state  vector  is  fundamental  to  state 
estimation  since  if  a  particular  state  is  unobservable,  its  value 
cannot  be  derived  from  the  output.  Alternatively,  a  combination  of 
states  may  be  observable,  thus  rendering  the  individual  states  indis¬ 
tinguishable  from  each  other.  Therefore,  for  many  feedback  control 
systems  it  is  assumed  that  the  entire  state  vector  is  available  through 
measurements.  Often,  however,  one  or  more  of  the  states  cannot  be 
observed  directly  with  only  the  measured  data,  and  an  approximate  state 
vector  must  be  substituted,  with  control  applied  to  the  estimate  instead. 

In  1963,  David  G.  Luenberger  derived  a  dynamic  device  which  would 
generate  estimates  of  the  inaccessible  plant  states  of  a  deterministic 
linear  system  from  exact  measurements  of  the  output  (Ref  14,  15). 

A  control  law  can  then  be  implemented  utilizing  the  resulting  full  state 
vector.  Again,  this  may  be  looked  upon  as  a  pole  placement  technique, 
and  as  long  as  the  noise  component  is  tiny  the  only  restriction  on  the 
arbitrarily-placed  observer  poles  is  that  they  have  faster  response 
characteristics  than  those  of  the  plant  matrix,  e.g.,  further  to  the 
left  in  the  s-plane,  for  rapid  convergence  to  the  true  state. 

For  the  system  under  consideration,  it  is  extremely  difficult  to 
determine  a  satellite's  position  in  the  orbit  through  conventional 
triangularization .  Only  line-of-sight  data  for  range  and  range- rate 
can  be  determined  to  any  great  degree  of  accuracy  due  to  the  alignment  of 
the  Earth,  Moon,  and  L3  point  and  the  size  of  the  orbit  (see  Figure  4) . 
However,  from  transponder  and  Doppler  shift  data,  accuracies  for  these 


radial  measurements  are  quite  impressive,  on  the  order  of  three  meters 
for  range  when  averaging  the  round  trip  time  of  radio  signals  or  down 


to  a  few  centimeters  for  laser  pulses  (Ref  18')  .  The  standard  deviation 
for  range-rate  is  .02  meters  per  second  ("Ref  S :  66 .  This  defines  the 
noise  in  range  and  range-rate  measurements  to  be  on  the  order  of  lO-^ 
so  that  double  precision  arithmetic  would  be  necessary  to  accurately 
simulate  the  effects.  Truncation  errors  due  to  computerized  numerical 
integration  are  large  enough  to  completely  obliterate  any  noise  from 
the  observations  themselves.  Therefore,  this  mav  be  considered  a 
strictly  deterministic  svstem. 

With  the  distances  involved,  range  measurements  varv  little  with 
nonradial  displacements.  For  example,  if  the  satellite  is  at  its 
maximum  y  position,  there  is  a  mere  .7%  variation  in  the  range.  This 
renders  the  plant  states  for  the  y  and  z  directions  essentially 
unobservable  and  estimates  must  be  constructed  using  a  Luenberger 
Observer. 

Repeating  the  state  equations  and  presenting  the  output  equations 
for  the  first  time  results  in  the  complete  system  to  be  considered: 
x(t)  =  A(t)x(t)  +  B(t)u(t) 

y(t)  =  C(t)x(t)  (36) 

The  output  or  measurement  vector,  v(t),  is  shown  in  component 
variational  form  in  (37)  .  The  general  range  and  range-rate  equations 


Range:  r(t)  =  J x(t)2  +  y(t)2  +  z(t)‘ 

.  -,  .  .  r(t) 

Range  Rate:  r(t)  =  v(t)  |?(t)T 


A 

^5x(t)2  +  6v(t)2  +  5z(t)^ 

[fix(t)5x(t)  +  6y(t)5y(t)  +  Sz(  t)  *  z  ( t)J  /  r(t) 


Figure  4.  Otate  Observability 


This  generates  an  output  matrix,  C(t),  of  the  following  form: 


36r(t) 

35r(t) 

36r ( t) 

36r(t) 

36r(t) 

3fir(t) 

3<5x(  t) 

3<SP  (t) 

35y(t) 

36P  (t) 
y 

36z(t) 

36P  (t) 

X 

3<5r(t) 

36r(t) 

36r(t) 

36r(t) 

36r(t) 

35r(t) 

35x(t) 

36P  (t)36y(t) 

36P  (t)36z(t) 

y 

36P  (t) 
2 

The  individual  C(t)  matrix  elements  were  verified  using  the 
definition  of  the  derivative.  For  example, 

3<Sr  _  [6r(i5x  =  2.001)  -  6r(6x  =  2)] 

11  36x  .001 


The  ranges  and  range-rates  were  examined  for  the  entire  orbit  and  were 
consistent  with  the  respective  data. 

Since  the  C(t)  matrix  is  a  2  x  6  matrix  rather  than  a  full  6x6 
matrix,  a  reduced  order  observer  will  need  to  be  considered.  The 
observer  state  and  output  equations  are 

xtt)  *  A(t)x(t)  +  B(t)u(t)  +  K(t)[y(t)  -  v(t)J 

y(t)  =  C(t)x(t)  (30) 

where  the  caret  represents  the  estimate  of  the  respective  vectors. 

Allowing  the  observer  reconstruction  error  to  be  defined  as 

e  (t)  =  x(t)  -  x(t)  (40) 

o 

gives 

e  (t)  =  A(t) x(t)  +  B(t)u(t)  -  A(t)x(t)  -  B(t)u(t) 
o 

-  K(t)[y(t)  -  y(t)]  (41) 

when  (36)  and  (39)  are  substituted  into  (40) . 


22 


Simplifying, 


eQ(t)  =  [A(t)  -  K(t)C(t)]e  (t)  (42) 

K(t)  is  termed  the  observer  gain  matrix,  the  dimensions  of  which  must 
be  6  x  2.  Elements  of  this  matrix  can  be  chosen  so  that  the  eigen¬ 
values  of 

j~A(t)  -  K(t)C(t)] 

will  all  lie  in  the  left-half  plane,  thus  allowing  the  steady  state 
value  of  e(t)  for  any  initial  conditions  to  approach  zero  as  t 
approaches  infinity. 

Since  only  one  of  the  three  system  modes  is  unstable,  modal 
control  theory  may  be  used  to  move  the  one  positive  pole  into  the 
other  plane,  following  the  same  derivational  steps  as  for  the  controller 

Through  the  modal  transformation  (27)  and  its  time  derivative 
(28)  the  system  equations  (36)  become 

n(t)  *  A*I(t)[A(t)A(t)-A(t)]n(t)  +  A-1  (t)B(t)u(t) 

y(t)  =  C(t) A(t)n(t)  (43) 

These  reduce  to 

n(t)  =  Jn(t)  +  A  1(t)B(t)u(t) 

y(t)  =  c(t)A(t)n(t)  (44) 

when  (19)  is  used. 

Substituting  the  equivalent  transformation  into  the  observer 
(39)  results  in 

A(t)*(t)  =[A(t)A(t)  -  A(t£)  n(t)  +  B(t)u(t) 

+  K(t)C(t)A(t)[n(t)  -  n(t)]  (45) 

which,  when  multiplied  by  A  1  ( t) ,  becomes 
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n(t)  =  Jn(t)  +  A  1(t)B(t)u(t) 

+  A_1(t)K(t)C(t)A(t)[n(t)  -  n(t)] 
Defining  the  modal  reconstruction  error  as 


e  (t)  =  n(t)  -  n(t) 
m 


gives 


e  (t)  =  n(t)  -  n(t)  (48) 

.  m 

Substituting  (43),  (45),  and  (47)  into  (48)  results  in  an  equation 
for  the  error. 

e  (t)  =  Jn(t)  +  A  !(t)B(t)u(t)  -  Jn(t)  -  A  1(t)B(t)u(t) 
m 

-  A-1(t)K(t)C(t)A(t)[n(t)  -  f|(t)] 
e  (t)  =  fj  -  A_1(t)K(t)C(t)A(t7]  e  (t)  (49) 

m  *—  m 

When  the  elements  of  the  observer  gain  matrix  are  chosen  so 


[j  -  A-1(t)K(t)C(t)A(t)J 


has  Poincare  exponents  with  only  negative  real  parts,  the  observer 
is  stable.  Thir-  will  give  identical  Poincare  exponents  when  trans¬ 
formed  back  to  the  physical  system. 

In  order  to  choose  the  elements  of  the  gain  matrix  easily,  the 
following  transformations  are  made: 

K"  ( t)  =  A-1(t)K(t) 

C'(t)  =  C(t)A(t)  (5C 


This  is  justified  since  the  observer  introduces  an  additional  set  of 
time  periodic  coefficients  and  Floquet  Theory  may  be  used  for  its 
analysis.  Now  (49)  becomes 

i  (t)  -  [j  -  K'(t)(T(t)]i  (t)  (51) 

m  _J  m 


t 


If  K'(t)  is  defined  as 


where  0  and  8  are  arbitrary  constants,  then  only  one  row  is  added  to 
J.  Thus,  by  appropriate  choices  of  0  and  8  the  third  root  can  be 
relocated  in  the  left-half  plane  as  desired.  Since  there  is  no 
analytical  method  to  conveniently  choose  gains  corresponding  to  desired 
pole  placement,  various  values,  both  constant  and  time-varying,  were 
implemented  and  the  resulting  Poincare  exponents  examined.  After  the 
appropriate  gain  has  been  determined,  the  K(t)  matrix  necessary  to 
stabilize  the  system  in  the  physical  variables  can  thus  be  easily 
derived  from  K^(t)  using  the  definition. 

Expanding  (51)  for  the  full  six-dimensional  error  vector  shows 


the  addition  of  only  one  row  to  influence  J. 


where 


J!  0  0  0 

0  0  0  0 

N  N  N  N 

2  3  4  5 

0  0X0 

4 

0  0  0  0 

0  0  0  u 

6 


Ni  -  —  1 1  ( t )  +  SC2i(t^ 

N2  =  -[flCizCt)  +  8C2  2  ( t)J 
n3  =  *3  ”  [j^c  1  3  ( t)  +  SC2  3  ( 
N4  -  -gcx4(t)  +  BC24(tT] 
Ns  =  - [0C 1 5 ( t)  +  8C25(tT] 
Ne  =  -|0Ci6(t)  +  8(^2  6  (tTj 


e  (t) 
ml 

e  (t) 
m2 

e  (t) 
m3 

i  e  (t) 
(  mu 


e  (t) 
m5 

e  (t) 
mg 


(53) 
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Now  an  explicit  equation  for  e  (t)  mav  be  formulated. 

m3 

V(t)  ■  '  08cn(t)  +  ec2i(t)3emi(ti  -|V12(0  +  BC22(o]e  (t) 

+  [\3  -  [jac- 3  (t)  +  eC23(t/]]em3(t)  -  pcl4(t)  +  SC2  4  ( tTj  ( t)  (54) 

-  [jlC^Ct)  +  BC^tjJe^t)  -  |7ci6(t)  +  BC26(t7]e  (t) 

Even  though  (54)  involves  the  entire  error  vector,  this  mav  be  simpli¬ 
fied  to 

%(t)  "  -  [^3(0  +  8C23(t3jem3(t) 

since  the  remaining  terms  exhibit  known  time  periodic  behavior  from 
(53)  and  so  may  be  combined  as  a  time-varving  particular  solution. 

Substituting  C 1 ( t)  for  Cj3(t)  and  £ 2 (t)  for  C23(t)  renders  the 
differential  equation  to  be  solved  as 

e  (t)  -  (x3  “  [«»Ci(0  +  652(t3}e  (t)  =  0  (55) 

Multiplying  both  sides  by  the  integrating  factor 
e-/[A3-[05i(t)  +  0C2(O]}dt 

gives 

e-/{x3-[0C1(t)  +  6£2(t)]]dt  ;  (t)  . 

m3 

C  -/fl 3  —  1  ( t )  +  6^2 ( 

[X3  -  [05!  (t)  +  B£2(tj]j  e  C  J  e  (t) 

The  solution  to  this  exact  differential  equation  is 
-/  fx3  -[0^(0  +  8£2 ( tm  dt 

e  (t)  e  J  *  c  (56) 

m3 

where  c  is  the  constant  of  integration. 


26 


Fourier  analysis  of  the  C^(t)  matrix  reveals  that  each  element 
is  made  up  of  a  constant  and  periodic  terms.  This  allows  £](t)  and 


C  2 ( c )  t0  be  expanded  into 

Ci (t)  =  Ci o  +  Cip(c) 

C2(t)  -  C20  +  C2p(c) 

where  the  o  signifies  constant  and  the  p  periodic  terms. 

Solving  for  e  (t)  and  substituting  for  Ci(t)  and  C-’(t)  gives 
m3 


e  (t)  =  ce 
m3 


/ Q  3  -  0CiO  -  ^Cip(t)  -  BC20  ~  BC2p(t)]  dt 
/Q3  -  0C 1  0  ~  SC2o]dt  -/[0Cip(t)  +  BC2p(tT]dt 


=  ce  e 

£3  -  0C 1  0  -  SC2q]t  -/[0Clp(t)  +  8C2o(  t)]  dt  (57) 

=  ce  e 

The  periodic  terms  oscillate  so  that  the  stability  characteristics 
are  dominated  by  the 


-  0Eio  -  SS20]  t 

e  ^ 

term.  Therefore,  the  Poincare  exponent  is  placed  by 

Xj  =  X3  -  0C 1 0  ~  3C2o 


(58) 


a  linear  form.  Here  X3  is  the  present  position  and  X3  is  the  desired 
position  of  the  Poincarfe  exponent. 
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V.  Controller /Ob server  Implementation 

The  purpose  of  designing  a  Luenberger  Observer  is  to  be  able 
to  feed  the  resulting  state  estimate  into  the  controller  which  would 
then  act  to  stabilize  the  svstem.  If  Che  observer  poles  are  to  the 
left  of  those  of  the  system,  keeping  in  mind  the  noise  component  re¬ 
strictions  metioned  earlier,  convergence  to  the  actual  state  vector 
is  rapid  and  the  controller  is  furnished  with  excellent  state  values. 
It  should  not,  however,  deteriorate  the  performance  of  the  controller. 

It  can  be  shown  that  observer  and  controller  eigenvalues  can  be 
assigned  independently,  thereby  creating  a  composite  and  vet  uncoupled 
set.  The  former  control  law,  (33),  modified  to  utilize  the  state 
estimate,  x(t),  becomes 

~1 

u(t)  =  -  F( t)x( t)  if  F(t)  =  -  GA3  (t) 

Inserting  this  into  (36)  results  in 

5(t)  =  A(t)x(t)  -  B(  t)  F(  t)  x(  t)  (59) 

y(t)  «  C(t)x(t) 

Recalling  the  reconstruction  error  in  phvsical  variables, 

e  (t)  =  x(t)  -  x(t)  (40) 

o 

and  substituting  into  (59)  above  yields 

x(t)  =  [A(t)  -  B(t)F(t)]x(t)  +  B(t)F(t)So(t)  (60) 

This  can  be  presented  in  partitioned  form  in  conjunction  with  the 
observer  reconstruction  error  equation  (42)  as 


x(t) 

A(t)  -  B(t)F(t) 

“1 

B(t)F(t) 

x  ( t) 

e  (t) 

L  0  J 

0 

A(t)  -  K(t)C(t) 

e  ( t) 
0 
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Now  defining  a  system  state  transition  matrix  as 


*s(t)  =  [*Jij]*s(o 

(62) 

where 

'V  =  A(t)  -  B(t)F(t) 

5  =  B(t)  F(t.) 

v  =  A(t)  -  K(t)C(t) 

Expanding 

the  partitioned  equations  yields 

$Sn(t>  =  i'<fSn(t)  +  ^S2l(t) 

(63a) 

$S12(t)  =  WSl2(t)  +  C<5S22(t) 

(63b) 

$S21 (t)  =  v$S2i (t) 

(63c) 

'fs22(t)  =  u<tS22(t) 

(63d) 

As  in  any 

state  transition  matrix,  the  initial  conditions  are 

$s(0)  =  I 

This  sets 

the  initial  conditions  for  the  submatrices  defined  bv 

the  partitions  as  well: 

$Sll(0)  =  I 

(64a) 

*si2(0)  =  [  0] 

(64b) 

$S2l(0)  =  [  0] 

(64c) 

$S22(0)  =  I 

(64d ) 

g 

(63c)  implies  that  4>  21  =  0  for  all  time,  so  (63a)  and  (63d)  become 


After  one  period  has  elapsed 


$S(t) 


$S] 2(t) 

I 

O 

$S22(t) 

(65) 


Forming  the  characteristic  equation  of  this  time-varving  composite 
system  gives 


fx  I  -  s 

I  plant  $ 


1 1  ( t ^  '^observer1'  -  = 


where  X  .  are  the  n  eigenvalues  of  the  svstem  and  -  , 

plant  •  obser' 

those  for  the  Luenberger  Observer.  This  separation  deior.'  . 
placing  an  observer  in  a  feedback  svstem  does  not  affect  it 
values.  The  two  sets  merely  produce  a  composite  formation. 
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VI .  Results,  Conclusions,  and  Recommendations 


Results.  Since  all  of  the  system's  states  cannot  be  measured 
directly,  a  suitable  approximation,  which  could  be  substituted  into 
the  control  law  developed  by  Captain  Shelton,  was  sought.  Therefore, 
a  logical  extension  to  his  work  was  the  implementation  of  a  Luenberger 
Observer  which  would  reconstruct  the  entire  state  vector.  As  developed 
in  Chapter  IV,  the  observer  moves  the  unstable  component  linearlv 
with  gain,  similar  to  the  case  of  the  controller.  Appropriate  elements 
of  the  modified  observer  gain  matrix  were  found  using  a  svstematic 
search  technique  and  are  depicted  in  Figure  5.  For  example,  a  value 
of  (3  =  2000.0  and  3  =  0.0  in  (58),  moves  the  root  from  2.392121  to 
-14.92828,  obviously  far  enough  into  the  left  half-plane  to  render  the 
system  stable.  The  slope  was  calculated  as  -0.00866. 

When  other  elements  of  the  modified  observer  gain  matrix  were 
filled  in  order  to  try  to  bring  the  oscillating  components  into 
negative  real  configuration,  the  linear  correlation  is  destroved. 

All  the  states  except  the  fourth  couple  and  it  becomes  virtual lv 
impossible  to  predict  the  effect  of  various  K'(t)  values.  Some 
examples  of  exponent-gain  correlation  are  presented  in  Table  ITT. 

It  was  eventually  possible  to  bring  the  oscillating  planar  mode  into 
negative  real  form  through  appropriate  application  of  gains,  but  to 
date,  the  same  has  not  been  accomplished  for  the  nonplanar  mode. 

The  best  has  been  to  move  the  real  part  slightly  farther  into  the 
left-half  plane,  e.g.  -0.045  ±  1.2321. 

With  no  nonplanar  motion  considered,  a  gain  matrix  producing  a 
set  of  Poincar4  exponents  which  fulfills  the  stabilization  criteria 


Table  III.  Observer  Exponent-Gain  Correlations 
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2.39  0.01 


without  experiencing  adverse  noise  effects  is 


0  0 

3.0  0 


1500.0  0  0 


The  corresponding  Poincare  exponents  are  -2.214,  -6.047,  -11.47, 

-2.39,  0.0  ±  1.1251,  with  -2.214  dominating  the  system's  behavior 
since  the  rest  are  further  in  the  left-half  plane. 

Examining  this  observer's  performance  by  evaluating  the  data 
retention  time  reveals  that  this  system  is  being  implemented  correctly. 

The  fastest  root,  -11.47,  exhibits  a  time  constant  of  approximately 
.1719  periods  or  5.07  days,  which  is  quite  sufficient  since  an  update 
is  taken  every  few  seconds.  The  model  under  consideration  is  strictly 
deterministic  with  observational  inaccuracies  inconsequential  when 
compared  to  those  imposed  by  numerical  integration.  It  is  even  difficult 
to  distinguish  noise  in  the  double  precision  mode.  Therefore,  this 
observer  is  not  in  any  danger  of  approaching  the  undesirable  behavior 
characteristics  of  a  differentiator,  and  is  the  basic  one  used  in 
various  test  cases. 

When  the  controller  and  observer  were  implemented  in  coni  unction, 
the  estimates  of  the  states  converged  rapidly  to  the  true  ones. 

Three  basic  cases  were  chosen  to  demonstrate  the  widest  range  of 
performance  characteristics.  The  first  (Case  A)  is  a  very  accurate 
observer,  assuming  a  14.9  m  accuracy  for  the  line-of-sight  direction  and 
a  14.9  km  accuracy  in  the  y,  with  orbit  initial  conditions  beyond  the 
range  of  the  controller's  ability  to  achieve  convergence.  The  second 
(Case  B)  was  chosen  as  very  close  to  the  periodic  orbit  so  that  controller 
behavior  could  be  checked  and  incorporated  an  observer  having  realistic 
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properties,  e.g.,  extremely  accurate  data  in  the  line-of-sight  direction 
and  almost  no  knowledge  of  cross  track  capability.  The  third  case 
(Case  C)  is  the  former  with  the  out-of-plane  mode  included  for  comparison. 

To  note  the  behavior  of  the  observer  alone,  the  case  of  no  control 
on  orbit  initial  conditions  (implying  exact  positioning  on  the  periodic 
orbit)  was  run  (Figure  6).  This  showed  the  observer  converging  fairly 
rapidly,  matching  the  true  states  exactly  by  at  least  the  end  of  the 
fifth  orbit. 

The  next  series  of  figures  (7-9)  depict  the  characteristic  time 
evaluation  of  the  modal  states.  With  the  extremely  accurate  observer, 
convergence  to  the  true  state  is  so  rapid  that  the  difference  between 
the  two  is  not  discernible  even  though  this  case  is  unstable.  The  other 
two  cases  depict  almost  total  convergence  by  the  end  of  the  third  orbit. 

In  addition,  appropriate  long-term  controller  behavior  is  demonstrated 
for  ni  and  r)3  in  that  these  elements  of  the  modal  state  vector  approach 
zero  with  time.  To  check  this  behavior  for  the  other  two  elements  would 
only  necessitate  a  longer  time  interval,  but  similar  tendencies  appear 
to  be  present.  The  z  mode,  when  affected  (See  Figure  9),  depicts  the 
expected  oscillatory  behavior. 

Figures  10,  11,  and  12  are  the  errors  between  the  modal  and  true 
states  with  time.  These  show  convergence  to  zero,  but  may  be  somewhat 
misleading  due  to  the  scale.  Approximate  exponential  envelopes  mav  be 
constructed  to  define  the  convergence  rate  more  closely  (Table  TV) . 

These  reveal  that  the  errors  are  dominated  by  the  largest  negative 
real  Poincar4  exponent  of  the  observer. 


Figure  6a. Uncontrolled,  bystem  with  Observer 


Controller  gain  0.0 

Orbit  initial  conditions 
x  -  0.0  km 
£x  ~  0.0  kn/sec 
y  =  0.0  km 
? y  =  0.0  kn/sec 
z  =  0.0  km 
=  0.0  kny^sec 

Five  orbits 


Observer  initial  conditions 
emi  =  Id. 96  m 
em2  =  14.96  m/sec 
em3  *  14959-97  km 
em4  =  14  959-97  km/s 

em5  =  0>0  kl\ 
en6  =  0.0  km/sec 

Observer  gain  elements 
Kh  =  3.0 


Examining  the  phase  portraits  for  tne  individual  modes  in 
each  test  case  (Figures  13,  14,  and  15)  reveals  tnat  once  again 
the  accurate  observer  of  Case  A  has  already  caused  the  estimate 
to  match  the  true  state  vector  identically  even  though  the  system 
is  inherently  unstable.  From  the  figures  it  is  somewhat  difficult 
to  determine  the  point  of  convergence  for  the  other  two  cases,  but 
the  listing  of  the  data  reveals  convergence  to  the  sixth  decimal 
by  the  third  orbit  in  Case  B  (nonplanar  mode  unchanged)  and  simi¬ 
larly  for  Case  C. 

Figure  14b.  is  especially  interesting  since  it  demonstrates 
the  suppression  of  an  unrealistically  large  initial  transient. 

The  observer  and  controller  have  widely  disparate  starting  positions 
and  for  the  first  four  points  (approximately  one-half  orbit)  the 
controller  actually  diverges  from  the  state  estimate.  Later,  the 
trend  returns  to  convergence  and  joint  motion  toward  zero. 

Figure  16  depicts  a  large  orbit  initial  condition  and  almost 
unrealistic  lack  of  knowledge  in  the  y-direction  observer  states. 
This  example  was  chosen  to  demonstrate  the  system’s  behavior  as  the 


Case  A; 


Case  B: 


Case  C: 


Table  IV. 


Exponential  Envelopes  for  Errors 
upper  envelope  =  1.833  x  lCT?  e 
lower  envelope  =  -1.183  x  10  '  e 

upper  envelope  =  1.843  x  10'^  e 

-4 

lower  envelope  =  -1.179  x  10  e 
upper  envelope  =  1.728  x  10  e 
lower  envelope  =  -1.185  x  10  e 


-2.214t 
-2.2l4t 
-2.214t 
-2. 214 t 
-2.214t 
-2.214t 
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observer  converges.  Even  with  such  a  massive  uncertainty,  the  ob¬ 
server  behaves  well.  Figure  17  was  included  to  exhibit  the  observer’s 
characteristic  of  converging  faster  with  increased  negative  posi¬ 
tioning  of  the  Poincare  exponent.  However,  the  data  reveals  that 
the  convergence  is  actually  slower  in  this  case.  This  occurs 
because  with  the  gain  combinations  utilized,  the  dominant  negative 
root  now  becomes  -1.01. 

Figures  18-21  describe  other  interesting  cases,  comparing 
effects  of  variations  in  observer  and  orbit  initial  conditions. 

Even  with  generally  large  initial  conditions,  trends  toward  state 
convergence  are  evident  after  only  one  orbit  (Figure  19)  and  more  so 
with  an  order  of  magnitude  difference  in  the  third  and  fourth  obser¬ 
ver  error  states  (Figure  20).  However,  in  this  instance,  as  with  an 
order  of  magnitude  improvement  in  its  initial  conditions  (Figure  21), 
the  controller  no  longer  appears  effective  unless  numerous  orbits 
are  necessary  for  convergence  to  occur. 

Conclusions 

These  various  cases  exhibit  all  the  points  this  study  was  meant 
to  investigate:  the  existence  of  a  stable  Luenberger  Observer  for 
the  restricted  four- body  problem  in  question,  the  observer  pole- 
placement  behavior,  which  proved  to  be  linear  as  in  the  case  of 
Captain  Shelton's  controller,  the  rapidity  of  convergence  to  the 
true  states,  and  the  long-term  observer  and  controller  behavior 
which  does  not  impede  controller  performance. 
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.Recommendations 


The  number  of  cases  which  could  be  tested  for  system  behavior 
evaluation  is  practically  unlimited,  and  almost  all  would  be 
valuable  if  resources  and  time  would  allow.  One  should  be  able 
to  determine  a  combination  of  gains  which  result  in  a  totally 
negative  real  system,  including  the  nonplanar  mode,  which  would 
then  be  evaluated  for  comparative  performance.  Hopefully,  the  long¬ 
term  oscillatory  behavior  could  then  be  reducea  or  eliminated. 

One  should  also  seek  analytical  alternatives  to  arbitrary  pole 
placement.  It  has  been  suggested  that  an  examination  of  3utterworth 
patterns  (Ref  12)  would  be  a  possible  route. 

An  extensive  evaluation  of  noise  effects  on  observer  and 
observer/controller  behavior  is  appropriate  as  well.  j£ven  though, 
operationally,  this  is  a  deterministic  system  almost  impervious  to 
noise,  it  is  essential  to  establish  the  left-half  plane  boundary 
for  the  positioning  of  the  observer  poles  in  order  to  avoid  unde¬ 
sirable  differentiator  properties. 

Another  interesting  characteristic  which  may  reveal  some 
useful  system  properties  is  the  presence  of  the  large  frequency 
component,  6.18,  evident  in  numerous  trial  and  error  gain  calculations. 
These  may  occur  with  either  positive  or  negative  real  parts  of  dis¬ 
similar  values  and  the  same  imaginary  part,  instead  of  the  expected 
complex  conjugate  pairs.  Perhaps  they  indicate  the  maximum  extent 
of  a  stable  observer  performance  region,  a  boundary  to  be  avoided, 
or  perhaps  the  damping  characteristics  of  such  a  root  would  render 
such  a  situation  favorable. 


Figure  7 a.  Modal  Vector  vs.  Time:  Case  A 


Controller  gain  0.5 
Orbit  initial  conditions 
x  =  4488  km 
px  -  4488  km /sec 
y  x  4488  km 
Py  =  4488  km/sec 
2  =  0.0  km 
2  =  0.0  km  /sec 


Observer  initial  conditions 
emi  =  14.96  m 
em2  *  14.96  m/sec 
em3  3  l4‘96  km 
en4  =  14.96  km/sec 
em5  =  0,0  km 
em6  =  0.0  km/sec 

Observer  initial  gain 

K ' ^2  =  3-6 


Five  orbits 


ETA  2  (o)  AND 

EST  JMRTE  (*)  *10’* 

-0.22  -0.1«  -0.06  0.02  0.10  0.16 


Figure  ?b.  Second  State 


44 


* 


00 


o.so 


t.00 


l.SO 

TIME 


2.00 


Figure  7c.  Third  State 
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Figure  8a.  Modal  Vector  vs.  Time:  Case  B 


Controller  gain  0.5 

Orbit  initial  conditions 
x  =  14.96  km 
Px  =  14.96  kra/sec 
y  =  14.96  km 
Py  =  14.96  km/sec 
z  =  0. 0  km 
Pz  =  0.0  km/sec 


Observer  initial  conditions 
@ml  =  14.96  m 
em2  =  14.96  m/sec 
em3  =  19959-9?  km 
em4  =  14959.97  km/sec 
em5  =  0.0  km 
em6  =  0<°  km/sec 


Observer  gain  elements 
K '  1 2  =  3.0 
K'3!  =  1500.0 


Five  orbits 


Figure  8c.  Third.  State 


E ST  III 


KTA6  (»>  AND 

EST IMRTE  (*)  *10“ 

-0 . 16  -0.12  -0.08  -0.04  0.00  0.04 


Figure  8f,  Sixth  State 


&TA1  (°>  AND 
EST  I  MATE  (»)  *10"* 

.12  -0.06  -0.04  0.00  0.04 


Figure  9a.  Modal  Vector  vs.  Time:  Case  C 


Controller  gain  0.  5 

Orbit  initial  conditions 
x  =  14.96  km 
Fx  =  14.96  km/sec 
y  =  14.96  km 
Py  =  14.96  km/sec 
z  =  14.96  km 
Pz  =  14.96  km/sec 


Five  orbits 


Observer  initial  conditions 
eml  =  lf-96  m 
em2  =  1^*96  m/sec 
em3  =  14959.97  km 

em4  =  14959.97  km/sec 
em5  *  14959-97  km 
em6  =  14959-97  km/sec 


Observer  gain  elements 
K'i2  =  3.0 
K'31  =  1500.0 


E ST  111 


Figure  10.  Error  vs.  Time:  Case  A 


Controller  gain  0.5 
Orbit  initial  conditions 
x  =  4488  km 
?x  =  4488  km/sec 
y  =  4488  km 
P  =  4488  km/sec 
z  =  0.0  km 
Pz  =  0,0  km/sec 


Two  orbits 


Observer  initial  conditions 
emi  =  1^.96  m 
era2  =  1^*96  m/sec 
em3  =  14.96  km 
em4  =  14,96  km/sec 
em5  =  0.0  km 
em6  =  0*0  km/sec 


Observer  gain  elements 
K* 12  =  3,0 
K  *31  =  1500.0 
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Figure  11.  Error  vs.  Times  Case  B 


Controller  gain  0.5 
Orbit  initial  conditions 
x  =  14. 96  km 
Px  =  14.96  km/sec 
y  =  14.96  km 
Py  =  14.96  km/sec 
z  =  0.0  km 
Pz  =  0.0  km/sec 


Observer  initial  conditions 
eml  =  14.96  m 
em2  =  14.96  m/sec 
em3  =  149*9.97  km 
en4  =  14959*97  km/sec 
em5  *  0.0  km 
em6  =  0.0  km/sec 


Two  orbits 
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Observer  gain  elements 
K’l2  =  3*0 
K*3i  =  1500.0 


ERROR 


Figure  12.  Error  vs.  Time;  Case  C 


Controller  gain  0.5 

Orbit  initial  conditions 
x  =  14.96  km 
?x  =  14.96  km/sec 
y  =  14.96  km 
ry  =  1^.96  km/sec 
z  =  14.96  km 
Pz  =  14  96  km/ sec 


Two  orbits 


Observer  initial  conditions 
eml  =  14.96  m 
em2  =  1^-96  m/sec 
em3  =  1^939 .97  km 
em4  =  14959*97  km/sec 
em5  =  14959-97  km 
em6  =  1^959-97  km/sec 


Observer  gain  elements 
K'l2  =  3.0 

K’3!  =  1500.0 
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Figure  13a.  Phase  Portraits:  Case  A 


Controller  gain  0.5 

Orbit  initial  conditions 
x  =  4488  fan 
Px  =  4488  fan/sec 
y  =  4488  fan 
Py  =  4488  fan/sec 
z  =  0.0  km 
Pz  =  0.0  km/sec 


Observer  initial  conditions 
eml=  14.96  m 
em2=  14.96  m/sec 
«m3=  14.96  km 
em4=  14.96  km /sec 
em5=  0.0  km 
em6=  0*0  km /sec 


Five  orbits 


Observer  gain  elements 


m2 


Figure  14a.  Phase  Portraits:  Case  B 


Controller  gain  0.5 

Orbit  initial  conditions 
x  =  14.96  km 
Px  =  14.96  km/sec 
y  =  14.96  km 
Py  =  14.96  km /sec 
z  =  0.0  km 
Pz  =  0.0  km/sec 


Five  orbits 


Observer  initial  conditions 
eml  =  14'96  m 
em2  =  14*96  m/sec 
em3  =  14959*97  km 
e,^  =  14959*97  km/sec 
em3  =  0.0  km 
era6  =  0.0  km/sec 


Observer  gain  elements 


k'31  =  1500.0 
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Figure  14c.  Fifth  and  Sixth  States 


69 


T 


3  *2013 

(D)  ZVlfl 


Figure  15a.  Phase  Portraits:  Case  G 


Controller  gain  0.5 

Crbit  initial  conditions 
x  =  14,96  km 
?x  =  14.96  km  /sec 
y  =  14,96  km 
Py  =  14.96  km /sec 
z  =  14.96  km 
P2  =  14.96  km/sec 


Five  orbits 


Observer  initial  conditions 

emi  =  14.96  m 

em2  =  14.96  m/sec 
em3  =  14959.97  km 
em4  =  14959.97  km/'sec 
em5  =  14959.97  km 
em6  =  14959.97  km/sec 

Observer  gain  elements 
K*12  =  3.0 

K*-Jl  =  1500.0 
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Figure  l6a.  Effect  of  Large  Observer  Initial  Conditions 


Controller  gain  0.5 
Orbit  initial  conditions 
x  =  4488  km 
Fx  =  4488  km/ sec 
y  =  4488  km 
Fy  -  4488  km/sec 
z  =  0.0  km 
r;z  =  0.0  km/sec 


Observer  initial  conditions 
eml  =  14- .  96  m 

em2  =  14.96  ra/seg 

em3  =  1*^96  x  10?  km 

em4-  =  ^  ,i+96  x  10°  km/ sec 

e  ^  =  0.0  km 
m3 

em6  =  6*6  km/sec 


'.  •/<*  or  cite 


Observer  gain  elements 
K*12  *  3.0 

K’ji  =  1500.0 
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Figure  l6c.  Third  State 


Figure  l6d.  Fourth  State 
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Figure  17a.  Case  of  Figure  16.  with 
Alternate  Observer  Gain  Matrix 


Controller  gain  0.5 
Orbit  initial  conditions 
x  =  4488  km 
Px  =  4488  km/sec 
y  =  4488  km 
Py  =  4488  km/sec 
z  =  0.0  km 
Pz  =  0.0  km/sec 


Observer  initial  conditions 

eml  =  14 -96  m 
em2  =14.96  m/sec 
em3  =  1 .4-96  x  10?  km 
em4  =  1.496  x  10°  km/se 
em5  =  0.0  km 
em6  -  0-0  km/sec 


Five  orbits 


Observer  gain  elements 
K'12  =  3.0 

K'31  =  4000.0 
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Figure  17c.  Third  State 
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Figure  17d.  Fourth  State 
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Figure  18a.  Additional  Example  of 
a  Faster  Observer 


Controller  gain  0.5 
Orbit  Initial  conditions 
x  =-50.42  km 
Px  =-50.42  km/sec 
y  =-50.42  km 
Py  =-50.42  km/sec 
z  =  0.0  km 
Pz  =  0.0  km/sec 
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Observer  initial  conditions 
eml  =  14.96  m 
em2  =  14.96  m/sec 
em3  =  1496.0  km 
em4  =  1496.0  km/sec 
em5  =  0.0  km 
em6  =  km/ sec 

Observer  gain  elements 
K'l2  =  3.0 
K#3!  =  2500.0 


One  orbit 
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Figure  19a.  Behavior  With  Large 
Initial  Conditions 


Controller  gain  0.5 
Orbit  initial  conditions 
x  =  4488  km 
Px  =  4488  km/ sec 
y  =  4488  km 
Py  =  4488  km/sec 

z  =  0.0  km 

P„  =  0.0  km/sec 
z 

One  orbit 

6:5 


Observer  initial  conditions 
eml  =  1-5  km 
em2  3  i.5  kry  sec 
®ra3  =  7.48  x  10-?  km 
em4  =  7.48  x  10*  km/sec 
em^  =  0.0  km 
em6  =  °.°  kra/sec 


Observer  gain  elements 
K ' 12  =  3.0 
K'3l  =  1500.0 
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Figure  20a.  Case  of  Figure  19.  With  Order  of  Magnitude 
Difference  in  Observer  Initial  Conditions 

Observer  initial  conditions 
eml  =  1-5  kn 
era 2  ~  1 . 5  fan /sec 
em3  =  7.48  x  10^  km 
era4  =  7.48  x  10°  km /sec 
em5  =  0.0  km 
em£  =  0.0  km /sec 
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